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ABSTRACT 

We discovered two simultaneous kHz quasi-periodic oscillations (QPOs) in the bright 
low-mass X-ray binary and Z source GX 5—1 with the Rossi X-ray Timing Explorer. In 
the X-ray color-color and hardness-intensity diagram a clear Z track is traced out, which 
shifted between observations. The frequencies of the two kHz QPOs increased from 
~215 Hz and ~500 Hz on the left part of the horizontal branch to ~700 Hz and ~ 890 
Hz, respectively, on the upper part of the normal branch. With increasing frequency the 
FWHM and rms amplitude (8.6-60 keV) of the higher-frequency kHz QPO decreased 
from 300 to 30 Hz , and from 6.6 to 2.4%, respectively. The FWHM and amplitude 
of the lower-frequency kHz QPO (50-100 Hz and 3-4% rms) did not correlate with the 
position of the source on the Z track. The kHz QPO separation was consistent with 
being constant at 298±11 Hz. Simultaneously with the kHz QPOs horizontal branch 
oscillations (HBOs) were detected with frequencies between 18 and 56 Hz. 

Subject headings: accretion, accretion disks — stars: individual (GX 5—1) - - stars: 
neutron — X-rays: stars 
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1. Introduction 

The bright low-mass X-ray binary (LMXB) and Z 
source (Hasinger & van der Klis 1989) GX 5—1 traces 
out a clear Z track in the X-ray color-color diagram 
(CD). The limbs of this track are called the horizon- 
tal branch (HB), normal branch (NB), and flaring 
branch (FB). It is thought (e.g., Hasinger & van der 
Klis 1989) that the mass accretion rate (M) increases 
when the source moves from the left end to the right 
end of the HB, down the NB, and onto the FB. Two 
types of quasi-periodic oscillations (QPOs) were so far 
detected in GX 5— 1 : on the HB oscillations called 
HBO were detected with a frequency between 13 and 
50 Hz (van der Klis et al. 1985), and on the NB oscil- 
lations called NBO with a frequency of ^6 Hz (e.g., 
Lewin et al. 1992). No QPOs were detected on the 
FB of GX5-1. 

Recent, two simultaneously kHz QPOs were dis- 
covered in the Z sources ScoX-1 (van der Klis 1996a; 

1997) , GX 17+2 (Wijnands et al. 1997b), Cyg X- 
2 (Wijnands et al. 1998), GX 340+0 (Jonker et al. 

1998) , and GX 349+2 (Zhang, Strohmayer, & Swank 
1998). The properties of these QPOs are remarkably 
similar to each other. They are best detected on the 
HB and the upper part of the NB. So far, Sco X-l 
is the only Z sources in which kHz QPOs have been 
detected on the lower part of the NB and the FB. 
The frequencies of the kHz QPOs increase from the 
left end of the HB to the upper part of the NB. In 
Sco X-l, their frequencies continue to increase as the 
source moves down the NB and onto the lower part of 
the FB. Similar kHz QPOs have been detected (e.g. 
Strohmayer et al. 1996; Berger et al. 1996; Zhang et 
al. 1996) in the less luminous LMXBs, called the atoll 
sources (Hasinger & van der Klis 1989). The many 
similarities (e.g., the frequencies and the dependence 
of the kHz QPO parameters on inferred M) between 
the kHz QPOs detected in the Z and atoll sources 
suggest that these QPOs are produced by the same 
physical mechanism in both types of LMXBs (see van 
der Klis 1998 for a review of kHz QPOs). 

In this paper we report the discovery of two si- 
multaneous kHz QPOs in the Z source GX 5—1. A 
preliminary announcement of the discovery of the kHz 
QPO was already made by van der Klis et al. (1996b). 

2. Observations and Analysis 

We observed GX 5—1 on 1996 November 2, 6, and 
16 (the AOl data), and on 1997 May 30, June 5, 



and July 25 (the A02 data) using the proportional 
counter array onboard the Rossi X-ray Timing Ex- 
plorer (RXTE). We obtained a total of 102 ksec of 
good data. During 17 % of the observing time only 3 
or 4 out of the 5 detectors were active. In constructing 
the CD and the hardness-intensity diagrams (HIDs) 
we used only the data of the three detectors that were 
always on. For the power spectral analysis we used 
all available data. Data were obtained in 129 photon 
energy channels (2-60 keV) with a time resolution of 
16 s. Simultaneous data were collected in two photon 
energy bands (2-8.6 and 8.6-60 keV; the AOl data) 
or four (2-5.0, 5.0-6.4, 6.4-8.6, and 8.6-60 keV; the 
A02 data) with a time resolution of 122 [is. The CD 
and HIDs were constructed using the 16 s data; power 
density spectra were calculated using the 122 /is data. 

To determine the properties of the kHz QPOs, we 
fitted the 64-4096 Hz power density spectra with a 
function containing one or two Lorentzians (to repre- 
sent the kHz QPOs), a constant plus a broad sinusoid 
(the dead-time modified Poisson level; Zhang et al. 
1995), and sometimes an extra Lorentzian with a cen- 
troid frequency below 50 Hz to represent the power 
spectral continuum below 200 Hz. The Very Large 
Event (van der Klis et al. 1997) window was set to 55 
/xs. Its effect on the Poisson noise was small and could 
be incorporated in the broad sinusoid. The proper- 
ties of the HBOs and the NBOs were determined by 
fitting the 0.125-192 Hz power density spectra with 
one or two Lorentzians (the HBO and its second har- 
monic; the NBO), and an extra Lorentzian (with a 
centroid frequency ^0 Hz) to represent the underly- 
ing continuum. 

The uncertainties in the fit parameters were de- 
termined using A% 2 = 1 and upper limits using 
A% 2 = 2.71, corresponding to 95% confidence lev- 
els. Upper limits on the kHz QPOs amplitude were 
determined using a FWHM of 150 Hz. When only 
one of the kHz QPOs was detected, the upper limit 
on the other kHz QPO was determined by fixing its 
frequency to the frequency of the detected kHz QPO 
plus or minus (depending on which kHz QPO was de- 
tected) the mean peak separation. Upper limits on 
the amplitude of the HBO and its second harmonic 
were determined using a FWHM of 10 Hz or 20 Hz, re- 
spectively, and the frequency of the second harmonic 
was fixed at twice the frequency of the HBO. 

In the CD and the HIDs, the soft color is defined 
as the logarithm of the 3.5-6.4 keV/2. 0-3. 5 keV count 
rate ratio, the hard color as the logarithm of the 9.7- 
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16.0keV/6.4-9.7keV count rate ratio, and the inten- 
sity is the logarithm of the 2.0-16.0 keV count rate. 
By using logarithmic values for the colors and the in- 
tensity, the S z parameterization (see Section ||) does 
not depend on the values of those quantities but only 
on their variations (Wijnands et al. 1997a). The 
count rates used in the diagrams are background- 
subtracted but not dead-time corrected. The dead- 
time correction is 3-5 %. 

3. Results 

The CD and HID of all data combined is shown in 
Figures |l|a and b, respectively. Clearly, in the HID 
(Fig. |l|b) the Z track is multiple. The reason for this 
is that it shifted slightly between observations. The 
AOl data form one Z track with an extra branch trail- 
ing off the FB (Fig. The A02 data are located 
at three different places in the HID. The 1997 July 25 
data (Fig. |l]e) show a HB and NB which are displaced 
from the AOl track to lower count rates. The upper 
NBs of the AOl and A02 data do not fall on top of 
each other, but the lower NBs do (Fig. |l|b), indicating 
that the NB/FB vertex moved around less than the 
HB/NB vertex. The 1997 May 30 and June 5 data 
cover a segment of the NB (Fig. |l]d) located slightly 
to the right of the 1997 July 25 NB. The rapid X-ray 
variability (see below) confirms that GX 5—1 was on 
the upper part of the NB during the 1997 May 30 
and June 5 observations. In the CD (Fig. |l|a) the Z 
track is clearly visible with the extra branch trailing 
the FB and a slight upturn of the HB, which is only 
marginally visible in the HIDs (Fig. |l|b) . The motion 
of the Z track is in the CD obscured by the errors in 
the X-ray colors. 

Because of the motion of the Z track in the HID, 
we decided to analyse the power spectra of the AOl, 
the May/ June A02 and the July A02 data (here- 
after referred to as the A02 data, unless otherwise 
mentioned) separately. We selected the power den- 
sity spectra according to the position of GX 5—1 on 
the Z track in the HIDs of the AOl (Fig. [j]c) and 
A02 (Fig. [ljs) data, and measured for each selection 
the average position of GX 5-1 along the Z track us- 
ing the S z parameterization (see e.g. Wijnands et al. 
1997a) applied to the combined CD (Fig. |l|a). The 
errors on the S z values are the standard deviation of 
each selection. 

Two simultaneous kHz QPOs were detected in the 
AOl and A02 HB and upper NB data (Fig. |). These 



kHz QPOs were most significant in the highest pho- 
ton energy band (8.6-60 keV). As adding the lower 
energy data did not improve the significance, we used 
the 8.6-60 keV energy band throughout our analysis. 
The rms amplitude upper limits and the detected val- 
ues of the kHz QPOs at different photon energies are 
displayed in Table |l|. Note, that not all kHz QPO 
detections were highly significant. For example, the 
three detections of the kHz QPOs in the AOl data 
with the highest S z value (see Fig. ||a) are only sig- 
nificant at the 1.7-2.3 a level. 

The properties of the kHz QPOs versus S z are 
shown in Figure |^. The two simultaneous kHz QPO 
peaks increased in frequency (Fig. ^a) from the left 
end of the HB (~500 Hz and ^215 Hz, respectively) 
to the HB/NB vertex and the upper part of the NB 
(~890 Hz and ~700 Hz). Over this interval the rms 
amplitude (Fig. ||c) and the FWHM of the higher- 
frequency peak decreased from 6.6% to 2.4%, and 
from 300 Hz to 30 Hz, respectively. The rms ampli- 
tude (3-4 %; Fig. |e) and FWHM (50-100 Hz) of the 
lower-frequency kHz QPO show no clear relation with 
S z . In the combined A02 1997 May 30 and June 5 
data kHz QPOs were detected only on the upper part 
of the NB, with an rms amplitude, FWHM, and fre- 
quency of 2.0±0.4%, 93t^Hz, and 557^ 4 21 Hz for the 
lower-frequency kHz QPO, and 2.4±g;|%, 122^Hz, 
and 856^7Hz for the higher- frequency kHz QPO, re- 
spectively. Lower down the NB, the upper limits on 
the kHz QPOs were 2.2% rms and 3.4% rms for the 
lower- and higher-frequency kHz QPOs, respectively. 
During the AOl, the A02 May/June and the A02 
July data the peak separation was consistent with be- 
ing constant at 297±17 Hz, 298±11 Hz, and 299±15 
Hz, respectively. The average peak separation for all 
observations combined was 298±11 Hz. 

Simultaneously with the kHz QPOs, the HBO and 
its second harmonic were detected. The frequency of 
the HBO increased from 18 to 56 Hz as the source 
moved from the left end of the HB to the HB/NB 
vertex and onto the upper part of the NB (Fig. ||b). 
Over this interval the rms amplitude of the HBO fun- 
damental decreased from 11% to 3% (8.6-60 keV; Fig. 

. The rms amplitude of the second harmonic (de- 
tected only below S z —1) decreased from 8% to 3% 
(Fig. |f). The FWHM of the fundamental increased 
from 5 to 24 Hz, and the FWHM of the second har- 
monic remained approximately constant at ~25 Hz. 
In the A02 May/ June data the HBO was also de- 
tected simultaneously with the kHz QPOs, with an 
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amplitude, FWHM, and frequency of 4.70±0.9% rms, 
16.0±0.8 Hz, and 50.9±0.2 Hz, respectively. NBOs 
with frequencies of ^5.7 Hz were detected between 
5 Z ~1.1 and S z ~1.8. The NBOs were detected simul- 
taneously with the HBO fundamental on the upper 
NB, but never simultaneously with the HBO second 
harmonic. No QPOs were detected further down the 
NB, on the FB, and on the extra branch. 

4. Discussion 

We have discovered two simultaneous kHz QPOs 
on the HB of the Z source GX5-1. KHz QPOs 
have now been discovered in all six of the known 
Z sources. The frequencies of both peaks increased 
with inferred M. The rms amplitude and the FWHM 
of the higher-frequency kHz QPO decreased as its 
frequency increased; the rms and the FWHM of 
the lower-frequency kHz QPO remained roughly con- 
stant. In Sco X-l the peak separation decreased when 
the source moved from the upper NB to the lower NB 
onto the FB. The peak separation for GX 5—1 is con- 
sistent with being constant, but also with being simi- 
lar to Sco X-l. This is true for the other Z sources as 
well, so that in the simplest description all peak sepa- 
rations vary, providing a challenge for beat frequency 
models where the peak separation is identified with 
the neutron star spin rate (see Psaltis et al. 1998 for 
a critical discussion of the peak separations). 

In the leftmost part of the HB, i.e. at the low- 
est inferred M, during the A02 observation we de- 
tected the lower-frequency kHz QPO at a frequency 
of 213.6l7 1 5 1 Hz. The corresponding higher-frequency 
kHz QPO has also a rather low frequency (504.5^ 12 ; 6 
Hz), but not unusually low (see e.g. Strohmayer et 
al. 1996; Ford et al. 1997; Jonker et al. 1998). 
The lower-frequency kHz QPO covers a total range 
in frequency of a factor 3.5. This QPO has the lowest 
frequency and the largest dynamic range for a kHz 
QPO so far detected in any LMXB. The frequency 
range falls in the 200-1000 Hz frequency range ex- 
pected if the lower-frequency kHz QPO is identified 
with the beat of the Keplerian frequency at the sonic- 
point with the neutron star spin frequency (Miller, 
Lamb, & Psaltis 1998). The corresponding range for 
the sonic-point radius would be 17-26 km, assuming 
a neutron star mass of 1.4 M Q and that the peak 
separation is the spin frequency. 

Recently, Stella & Vietri (1998) proposed that the 
low frequency QPOs observed in atoll sources are due 



to a precession of the innermost disk region, domi- 
nated by the Lense-Thirring effect. They also sug- 
gested that the HBO in Z sources could arise from 
the same mechanism. However, the precession fre- 
quencies predicted for GX 5—1 are about a factor 2.5 
smaller than that of the HBO frequencies (see also 
Stella & Vietri 1998), unless I/M (with I the moment 
of inertia of the neutron star and M its mass) is a 
factor of 2.5 higher than predicted by neutron star 
models with realistic equations of state, the kHz QPO 
peak difference is half the neutron star spin frequency, 
or the HBOs are not the fundamental frequencies but 
the second harmonics. For a detailed discussion about 
the Lense-Thirring interpretation of the HBO in Z 
sources we refer to Psaltis et al. (1998) 

The Z track of GX 5—1 shifts between observa- 
tions. The Z track of the AOl data is at ~4.5% higher 
count rates and ~2% softer hard colors than the Z 
track of the A02 data. Similar shift were already 
reported by Kuulkers et al. (1994) using EXOSAT 
data, although in their observations the variations in 
count rate and colors were slightly larger (~8% and 
~6%, respectively). When plotting the frequency of 
the higher-frequency kHz QPO and the HBO versus 
count rate (Fig. ||a and c) the AOl data are slightly 
displaced to higher count rates compared to the A02 
data. If we correct the count rates of the AOl data 
for the shifts in count rate derived above the kHz 
QPOs and HBO frequencies as a function of count 
rate are consistent with being the same during both 
epochs (Fig. ||b and d). The same is also true for 
the frequency of the lower-frequency kHz QPO (not 
shown here) and the rms and the FWHM of all the 
QPOs. The only difference could be the rms ampli- 
tude of the higher-frequency kHz QPOs at the lowest 
inferred mass accretion rate, however, the difference 
is marginal. This demonstrates that not the X-ray 
count rate determines the properties of the QPOs 
(both the kHz QPOs and the HBO), but the posi- 
tion of GX 5—1 on the Z track does. Kuulkers et al. 
(1994) also reported that the HBO properties were 
consistent with being the same when the Z track of 
GX 5—1 was displaced in the HID. In their observa- 
tions the tracks were displaced more than our tracks 
indicating that even for larger shifts the HBO prop- 
erties do not change. Jonker et al. (1998) found the 
same for the HBO (the statistics did not allow to con- 
clude anything for the kHz QPOs) in GX 340+0. For 
Cygnus X-2 (Wijnands et al. 1997) it was shown that 
the NBO properties changed significantly when the Z 
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track of this source moved around in the HID. How- 
ever, due to limited available data nothing could be 
concluded for the HBO properties in Cygnus X-2. 
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Table 1 
Amplitude of the kHz QPOs a 



Energy Lower-frequency pcak b Higher-frequency peak c 
(keV) (rms) (rms) 



2-5.0 
5.0-6.4 
6.4-8.6 
8.6-60 


< 1.3% 

< 2.2% 

< 2.3% 
3.1±°;°% 


< 2.0% 

< 2.6% 
2.6±0.5% 
4.7±0.4% 


a The values 


were calculated at Sz, 


=0.86±0.05 


b FWHM = 


98I34 Hz, frequency = 


= 340±^ Hz 


C FWHM = 


150^24 Hz, frequency 


= 647±10 Hz 
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Fig. 1. — The color-color and hardness- intensity dia- 
grams of GX 5-1. (a) the CD of all data, (b) the HID 
of all data, (c) the HID of the AOl data, (d) the HID 
of the A02 1997 May 30 and June 5 data, and (e) the 
HID of the A02 1997 July 25 data. In all diagrams 
the soft color is the logarithm of the 3.5-6.4 keV/2.0- 
3.5 keV count rate ratio, the hard color the logarithm 
of the 9.7-16.0 keV/6.4-9.7 keV count rate ratio, and 
the horizontal axis the logarithm of the 2.0-16.0 keV 
count rate. Background was subtracted but no dead- 
time correction (3-5%) has been applied. All points 
are 16 second averages. Typical error bars are plotted 
in the bottom right corner of (a) and (b). The gray 
solid line in (a) is the track used to calculated the S z 
values (see Section ||). 
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Fig. 2. — Typical Leahy normalized power density 
spectrum in the energy range 8.6-60 keV showing the 
two simultaneous kHz QPOs. 
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Fig. 3. — The properties of the kHz QPOs versus pos- 
tion on the Z track, S z (see Section ||) ; (a) frequency 
of the kHz QPOs, (b) frequency of the HBO funda- 
mental and 2nd harmonic, (c) the rms amplitude of 
the higher- frequency kHz QPO, (d) rms amplitude 
of the HBO fundamental, (e) rms amplitude of the 
lower-frequency kHz QPO, and (f) rms amplitude of 
the HBO 2nd harmonic. The black points are the 
A01 data and the red points the A02 July 25 data. 
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Fig. 4. — The higher-frequency kHz QPO frequencies 
(a and 6) and the HBO frequencies (c and d) versus 
the intensity. The definition of the intensity is the 
same as m Fig. 0. The bullets are the AOl data and 
the crosses the A02 July 25 data. In b and d the 
count rates of the AOl data were corrected for the 
shift of the Z track (see Fig. Section ||) . 
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